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Constitutive equations for the long- wavelength behaviour of the orientational dynamics of a super- 
cooled liquid are derived using a projection-operator technique and resulting expressions for light- 
scattering spectra are formulated. We thus extend recent studies for axially symmetric molecules to 
the general case of arbitrarily shaped rigid molecules. The second part of the discussion considers 
hydrodynamic energy-fluctuations and thus arrives at expressions for light-scattering spectra which 
also include a Rayleigh-line. The role of the memory-kernels in the theory is treated in detail. In 
particular, the derivation of a theory that formally resembles earlier approaches to the problem 
is presented using a mathematically rigorous description of the Laplace-transforms of correlation- 
functions. 

PACS numbers: 64.70.Pf, 78.35+c, 61.25.Em 



I. INTRODUCTION 

Light-scattering methods have proven to be a power- 
ful tool in the experimental investigation of super-cooled 
liquids, especially glass-forming substances, see for in- 
stance the works by Dreyfus and Pick [1-3] or those by 
Cummins [20]. In particular, the cross- or depolarised 
spectra, i.e., measurements in the VH-conhguration, have 
revealed many details of the process of structural re- 
laxation while cooling such liquids. As a glass-forming 
system approaches the glass-transition, the cross-over 
from pure liquid-like to solid-like behaviour is displayed 
in the scattering-spectra. Mode-Coupling-Theory for 
the glass-transition has helped to understand many of 
the phenomenological features of these spectra and sup- 
ported the interpretation of the underlying physical pro- 
cesses. However, a complete theoretical formulation of 
the physics of light-scattering in super-cooled liquids 
has not been available until recently. There are theo- 
ries which properly describe the Rytov-Dip in the high- 
temperature regime, which is due to the relaxation of 
shear-modes, i.e., a coupling of translational and rota- 
tional dynamics. These results stress the importance of 
including the rotational dynamics of the liquid in the 
theoretical description. On the other hand using mode- 
coupling-theory one can explain the appearance of the 
shear-waves in the dynamics at lower temperatures in 
place of the diffusive Rytov-Dip. The cross-over regime, 
however, could not be described. Recently, several stud- 
ies have been published, where phenomenologically ob- 
tained generalised Langevin equations have successfully 
been used to explain the light-scattering spectra of super- 
cooled liquids [1-4] . A microscopical derivation of these 
equations using the long- wavelength expansion of an ap- 
plication of the projection-operator formalism of Mori 
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and Zwanzig has been given for some special cases in 
[5, 6]. Although the first attempt to fully describe the ro- 
tational dynamics using a projection-operator technique 
was already made by Anderson and Pecora in 1971 [7], in 
[6] it was shown that their theory is not suited to describe 
the cross-over regime properly. Anderson and Pecora did 
not allow for a frequency-dependence of the transport- 
coefficients in their theory, and their choice of variables 
is not suited for the application of generalised hydrody- 
namics [8], which is the proper framework to treat the 
problem and within which the present results have been 
derived. Our results also rely on the tensorial formalism 
developed by Schilling and collaborators [9-11] and gen- 
eral results for the form of the light-scattering spectra 
due to symmetries [12]. However, the available models 
have simplified the physical system in question by only 
considering axially symmetric molecules and neglecting 
fluctuations of the local energy density. It is the aim 
of the present work to give a complete description of 
the rotational dynamics in the long-wavelength regime 
by applying the projection-operator technique of Mori 
and Zwanzig to the problem. The focus is mainly laid 
on a formal and rigorous derivation and discussion of the 
equations. We can therefore extend the above results to 
the cases of arbitrarily shaped rigid molecules and finally 
also include energy-fluctuations in the theory. 

The paper is organised as follows. In section II we de- 
fine the distinguished variables of our theory and discuss 
their basic properties as far as they are relevant for the 
application of the projection-operator formalism. Section 
III is concerned with the derivation of the generalised 
Langevin equations of the system in the time and the 
complex-frequency domain. The equations will be solved 
in section IV, where also the light-scattering spectra are 
computed and discussed. Some remarks regarding the 
positivity of the spectra will be made in section V. An 
extension of the theory that includes fluctuations of the 
energy-density is given in the following section. To be 
in accordance with earlier approaches to the problem we 
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will finally show in section VII how a theory that only 
considers fluctuations of the dielectric tensor can be de- 
rived from ours. 



II. HYDRODYNAMIC VARIABLES 

We consider a classical, super-cooled liquid consisting 
of N rigid molecules of arbitrary shape. The molecules 
are optically inactive and there is no orientational order 
present in the system. The liquid is enclosed in a volume 
V, a closed compact and simply connected subset of R 3 
with suitably smooth boundary, and held in equilibrium 
at a temperature T. We anticipate the thermodynamic 
limit N, V — > oo with the average density n :— N/V 
being held constant. The molecules have translational 
and rotational degrees of freedom, which are centre-of- 
mass positions R a and Euler-angles Cl a respectively. The 
Eulcr-anglcs describe the relative orientation of the prin- 
cipal axes of the tensor of inertia or the polarisability 
tensor with respect to the laboratory frame of reference. 
The two interpretations only differ by a global rotation 
which the system is invariant under. The phase-space is 
identified with R 6JV + 6JV i rp ne dynamics f ^ ne S y S tem is 
driven by the its Poisson-brackets or Liouvillc-operator 
respectively, A = {H, A} =: iCA, where H denotes the 
Hamilton-function of the liquid, 
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The symbol P a denotes the ath particle's translational 
momentum, J a ^ labels the ith component of the angular 
momentum of particle a and Ii is the respective eigen- 
value of the tensor of inertia. The relation between the 
moments of inertia and the canonical rotational momen- 
tum is given in [13]. Thus the volume- element of phase- 
space reads 



dT = II dR a dP a dn a dJ 



We assume the interaction energy V only to depend on 
positions R a and Euler-angles fl a , but neither on trans- 
lational nor orientational momenta. As a phase-space 
variable wc understand a function A : T — > U where 
U is a finite dimensional real or complex vector-space 
such that A is an element of L*(T) n L 2 p (T). The sub- 
script p indicates integration with respect to the equi- 
librium measure exp(— f3H)Z~ 1 , i.e., thermal averaging 
(•). Phase-space variables may also depend on an ad- 
ditional parameter, say r e V, and are assumed to 
have a Fourier-transform in this parameter. We will 
only consider fluctuations of phase-space variables, 8 A := 
A — (A), where in the following the 5 is implied by an 
abuse of notation. Statistical correlations are written in 
terms of the Kubo-scalar-product (A | B) := (SA*5B). 
With this construction all integrals are well-defined and 



the space of phase-space variables becomes a Hilbert- 
space. Parameter-dependent variables will further be 
represented by their corresponding spatial Fourier-modes 
A(q) := J v e lq ' r A{r) dr. Since we are free to define the 
orientation of the laboratory frame we will choose q par- 
allel to the system's z-axis. The variables we assume to 
exhibit hydrodynamical behaviour and will therefore be 
treated using a projection-operator approach comprise 
the fluctuations of the mass-density 
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and the mass-current, 



J(q) :- 



1 N 



mnv(q), 



(11.1) 



(II.2) 
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and variables describing collective rotational motion. By 
definition, the Hamilton-function is invariant under spa- 
tial rotations which favours a description relying on ir- 
reducible representations of the rotation-group 50(3), 
i.e. so-called spherical coordinates. For each solution of 
the equations of motion its whole orbit under the group 
action also is a solution. Thus every solution carries a 
representation of the group. The finite irreducible repre- 
sentations of SO(3) are unitary linear maps V which are 
classified by rank j. The helicity v G [— j, j] n Z denotes 
the ^th component of each vector in such a subspace. 

In this notion the mass-density is a scalar variable 
(j = 0), and the mass-current a vector (j = 1). To 
describe the orientation of a molecule and the quantities 
involved in light-scattering, namely the particle's polar- 
isability tensor [14], we will use a j — 2 object, i.e., a 
symmetric traceless tensor of rank two. We consider a 
single molecule's polarisability tensor Q^, a j — 2 ob- 
ject with helicity m, which is given in a certain body- 
fixed coordinate-system, rotate it back into the labora- 
tory frame of reference and add up the contribution for 
the whole liquid. Thus we find as a variable describing 
the relevant orientational properties 

1 N 

:= 7^ E e iqR "VW{n a )Q», (II.3) 

V JV a= l 

where the quantities 25^(0) are irreducible represen- 
tations of the rotation-group of rank j = 2, rotation- 
matrices on M 5 . (NB: Wc adopt the summation conven- 
tion for all latin- typeset indices.) This form of the vari- 
able and the fact that all particles have the same local 
properties suggest the use of rotation-matrices T>) l m(Q, a ) 
as the primary objects of interest. The definition of those 
follows [13]. To be in accordance with other theoretical 
work in the field, e.g. [5, 10], we define the orientational 
density 

P,m(q) ■= i 2 ^Y, eiqRaV f^)- (H.4) 
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In the small g-limit the time-derivative of this density 
will also be considered hydrodynamic. This assumption 
is important for proper modelling of the memory-kernels 
as was already pointed out in [6] . We define 

jfim : ~ Pfim (II.5) 

and call it an orientational current, although it is not de- 
rived from a conservation law. Finally, we will also con- 
sider long- wavelength energy- fluctuations. But to keep 
consistent with [4, 6] we will neglect them for the time be- 
ing, which is justified due to the closeness of the Landau- 
Plazcek ratio to unity. The necessary extensions to the 
theory will be presented in section VI. A treatment of 
general tensors of arbitrary rank is also possible. Exten- 
sive formal work has been done by Schilling et al. [10, 15]. 

A. Static Averages 

Static correlation-functions, which determine orthog- 
onality and normalisation of the dynamical variables, 
can be evaluated using the rotational invariance of the 
Hamiltonian function. According to the selection rules 
presented in Appendix B the q — averages of two func- 
tions are diagonal with respect to both helicity and rank. 
For finite q one may expand the correlator in terms of 
the wave- vector about the value q = 0. According to 
the selection-rules of Appendix A fixing the wave- vector 
parallel to e z , i.e., in spherical tensor notation q = qeo 
achieves decoupling of phase-space variables with respect 
to helicity. Then the mass-density-orientation correlator 
has the first non- vanishing term at order q 2 , cf. (A.6). In 
the long-wavelength limit mass-density and orientational 
density may therefore be considered orthogonal. Also, 
we can incorporate time-inversion-parity and establish 
orthogonality between densities and currents. 

The correlator between the orientational currents and 
the mass-current (p^m | J M ) is identically zero. This is 
proven easily. The orientational current breaks up into a 
translational part, which is linear in _P, and a rotational 
part which contains the derivatives of the Euler- angles. 
The interaction energy in the Hamilton-function is in- 
dependent of the momenta. Therefore, the correlation 
of the translational part of j Mm with the mass-current 
involves taking the expectation-value of the Wigner- 
functions themselves which by the assumptions on the 
system vanishes. The correlation of the rotational part 
of the orientational current with the mass-current van- 
ishes due to gaussian averaging over the translational mo- 
menta in the mass-current. Hence in the long- wavelength 
limit the hydrodynamical variables mass-density, mass- 
current, orientational density, and orientational current 
are mutually orthogonal. 

Their normalisations read, evaluated at q = 0, 

Smn '■— (Pfim \ Pfm) ari d Vl mn '.— (p^ m \ p^n) ■ (H.6) 

In Appendix B we will show that these matrices are 
positive definite, and hence invertible. Using the static 



structure factor S{q —> 0) and the associated isothermal 
sound- velocity Sc 2 := —: v 2 we have 

7Yl 2 V 2 

(p\p) = — — and {Jii{Q)\Ju{q))=5 IJtv m 2 v 2 . 

B. Equations of Continuity 

Both the mass-density and its current obey equations 
of continuity. The conservation law for the mass-density 
in spherical coordinates is the same as in cartesian coor- 
dinates 

p(q, t) = iq_ m J m (q, t). (II. 7) 

The conservation of momentum is achieved by using the 
stress tensor n which is a tensor having only its j = and 
j = 2 components non-vanishing, because the cartesian 
representation of n is a symmetric matrix. The conserva- 
tion law can be derived from the cartesian version, cf. for 
instance [6], using the known transformation laws [13], 

j„(g) = ^(lM0|j^)n(%). (II.8) 

3 

The coefficients of the linear combination in (II. 8) are 
Clebsch-Gordan-coefficients. Also, the relation q = qeo 
has been used here. Therefore the sum contains only 
terms for j = 0, 2. The scalar part, i.e., the j = com- 
ponent, is interpreted as pressure, whereas the tensorial 
part, i.e., the j = 2 components, describe shear. 

III. THE PROJECTION-OPERATOR 
TECHNIQUE 

In this section, we will derive the equations of mo- 
tion for the variables introduced before by applying 
the projection-operator technique of Zwanzig and Mori 
[16, 17]. The thus obtained equations have the form of 
generalised Langevin equations. The coupling between 
the conserved quantities and the orientational degrees of 
freedom will be achieved through the use of frequency- 
dependent memory-kernels which will have the function 
of transport coefficients. A light-scattering experiment 
only probes the dynamics of the liquid on a spatial scale 
which is large compared to any microscopic length-scale. 
It is therefore sufficient to derive the equations in leading 
order of the wave-vector q. Although the spatial scales 
are well separated the time-scales are not. It will turn 
out that by using frequency dependent memory-kernels 
the theory is able to describe the dynamics on a wide 
range of scales they as can be observed experimentally. 

A. Projection Operator 

In order to close the system of equations (II. 5, II. 7, 
II. 8) we need to derive constitutive equations for n(°), 
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II^ 2 ), and the orientational force p nn . As pointed out be- 
fore we will apply a projection-operator technique. Then 
the dynamics will naturally be split into three parts, the 
intrinsic dynamics of the hydrodynamic subspace, the 
mixing dynamics and those in the orthogonal comple- 
ment. The latter will generally be interpreted as the 
'bath degrees of freedom', because their correlation to the 
distinguished variables vanish. The variables onto which 
we will project the dynamics have already been intro- 
duced in section II, leading to the following projection- 
operator 

V(q) :=^(Omn \p»m{q)){p»n{q)\ 



(III.l) 



+ \p(q))(p(q)\; 



E 



\Mq))(M<i)\ 



We note that V is a projection only up to order 0(q 2 ), 
which is the convenient setting for the long-wavelength 
studies discussed in this paper. The time-evolution gen- 
erated by the Liouville-operator is given by an opera- 
tor R(t) which can be rewritten using the projection- 
operators V and Q := 1 — V, [6] 

R(t)=R(t)V + f R(s)ViCR'(t-s)ds + R'(t). (III.2) 
Jo 

A reduced time-evolution-operator has thereby been de- 
fined by 

R'(t) := Qe lQCQt Q. 

Equation (III. 2) is completely equivalent to the gener- 
alised Langevin equations discussed in [16, 17]. For each 
phase-space-variable A the evolved variable R'A is or- 
thogonal to each of the hydrodynamic variables. This 
property identifies it as noise and justifies its neglect for 
the remainder of this paper. 

We can simplify the integral-kernel ViLR! by operating 
with C to the left, 



V{-iC)R' = Y.i^Un |/W<Z)> (hn{q)\R' 



E^I J M('?))(n(f' ) (9)l(iM,io|jM)^- 



(III.3) 



The application of the evolution-operator (III. 2) to stress 
tensor and orientational force yields the missing equa- 
tions of motion. However, the price to be paid is the in- 
troduction of frequency dependent memory-kernels. Be- 
cause the reduced dynamics do not show any hydrody- 
namic poles, these kernels can be approximated by their 
q = value. 



A light-scattering experiment measures the Laplace- 
transform of the autocorrelation of the dielectric tensor 
or more precisely the rcal-frcqucncy-limit of its imagi- 
nary part [14]. A formulation of the equations of mo- 
tion to directly yield these correlator is more useful than 
the formulation in terms of generalised Langevin equa- 
tions. Also, the integro-differential equations will reduce 
to matrix-equations which can easily be solved. However, 
the equations in the time-domain cannot be derived from 
the matrix counterparts as both the Laplace-transform 
and canonical averaging are not invertible. We will there- 
fore perform all calculations in the complex-frequency 
domain, but also present the equations of motion in the 
time-domain, where the individual couplings between the 
various physical quantities can readily be read off. 



B. The Correlation Matrix 

According to [16] the Laplace-transform of a correla- 
tion matrix C(z) of a vector of hydrodynamic variables 
A = (Ai, . . . , A n ) obeys the following equation 



C(z)-S^ Sz-i{A\A) + {A\ R'A) ■S~ 1 = -Id, 

(III.4) 

where Id denotes the unit-matrix. All variables appear 
in their Laplace-transformed form which is convention- 



The symbol 
:= (A^Aj). 



ally set to £{f(t)}(z) := i j* e izt f(t) dt. 
S denotes the static correlation-matrix SV 
Given this formula we only need to evaluate a few static 
and dynamic correlation-functions to set up the equa- 
tions of motion. Due to time-inversion parities and rota- 
tional invariance of the Hamilton-function the only non- 
vanishing static correlators are 



{Pftm | Pan) = -^mn and (p \ J ) 



-mV. (III.5) 



The dynamical correlation-functions involving the re- 
duced time-evolution operator R'(t) are memory- kernels 
which can be related to transport-coefficients using the 
so-called Green-Kubo relations [16]. When calculating 
the light-scattering spectra we can identify the transport- 
coefficients by comparing with the theory of the simple 
fluid. Most of these coefficients are of order 0(q 2 ) and 
vanish in the long- wavelength limit. If we set Ai to one 
of those densities whose current is an element of the hy- 
drodynamic subspace the reduced dynamical correlator 
will vanish, because then QA = 0. This applies to all 
correlators involving p and p^m- Setting Ai = J a the 
correlator can be rephrased using the equation of conti- 
nuity (II. 8) and results in the autocorrelations of pressure 
p := -^in^ ) and the shears n^ 2 '. For the non-vanishing 

elements of the matrix (A | R'A) we find 

(p(q)\R'(z)p(q)) =:f,o{z) + 0(q 2 ) 

defining a bulk viscosity The correlation between 
orientational forces and pressure defines translation- 
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rotation couplings 

(p, n | i?'(z)n( 2 )) =: -S^p n (z) + 0(q 2 ). 
The autocorrelation of shear reads 

(U^( q )\R'(z)U^(q)) =:Uz)+0(q 2 ). 

Finally, the reduced autocorrelation of the orientational 
forces defines a rotational friction 

(Ptm(q) I R'(z)p^ r {q)) =: t nr (z). 

With these definitions and the decoupling of correlators 
with respect to helicities (A.l) we find three matrices rep- 
resenting the operator M(z) := Sz - i (A \ A) + (A \ R'A) , 
one for each modulus of helicity. For notational conve- 
niences we define some shorthands 

C v := (Iv, 10|2i/) and fj := C$fj a + fj b . (III.6) 

Using block-matrix notation and A := (p ur , p ur , J v , p) 
where r G \—v\v\C\7L we arrive at 



zS nr 



i/=±2 



(III.7) 



Z S flV ^ 71 T ^ 

M v =±i = I —iVL nr zQ nr + f nr -iqC±ip n 

;„n..r,T r>2 „2~ _l™2„,2. 



iqC±ip I r C\ x q 2 f] s + m 2 v 2 Zj 



(III.8) 



( zS nr 


V o 



iqC a pJ 






-iqC'oPr, 



q 2 f] + zm 2 v 2 qm 2 v 2 



(III.9) 

In this form, the equations of motion (III. 4) can be solved 
and the scattering spectra written down in a microscopi- 
cally exact form. There is an alternative approach going 
the route via constitutive equations in the time domain. 
The explicit procedure has been demonstrated eg. in [6]. 
Here, we will just give the results for our set of vari- 
ables without explicitly mentioning the noise term. A 
couple of abbreviations are introduced. First we absorb 
the frequency and structure-factor matrices into the cor- 
responding memory- kernels, 



and define 



Pm i T ) ■= (H 1 )mnPn(T), 
rmr(^) • )mnJ^nr (^")) 

(A')mr := (0) r 



knT 



(111.10) 
(III.ll) 

(111.12) 
(111.13) 



The viscosities with and without tilde are related to each 
other by a factor of jr^y • The equations for the mass- 
current turn out to be 

j±i{l,t) = iqC±i {p m * p±L >m (q)) (t) 

-£-C 2 ±1 (Vs*J±i(q))(t) (IH.14) 



and 

j (q,t) = iq 



c 2 p- 



iq 



iqC 



iqC'o 



(»m*tii(q)) (t) + l -^(r, s *Mq))(t) 



(111.15) 



where the asterisk denotes convolution in time. There- 
fore, using (II. 7), the equation for the mass-density is 

a 2 

p(q, t) = -q 2 c 2 p{q, t)-^-(n* p(q)) (t) 
mn 

- q 2 C (p m * P0,m 

(q))(s). (111.16) 

The equations for the orientational densities, with v = 
0, ±1 are 

p vr {q,t) = -Ld mr p» m (q,t) - (T mr * pvm{q)) (s) 

- (A'Wa— (Pm * Mq)) (a). (IH.17) 
mn 

And for v = ±2 one finds 

P±2,r(<7, t) = -U mr p±2, m {q, t) - (T mr * p±2, m (q)) (s). 

(111.18) 

These equations have the same functional form as the 
ones presented in [6] . Due to the use of Wigner- functions 
instead of a single orientational variable they form a 
vector-valued theory in contrast to the scalar theory de- 
veloped in [6] . For a comparison of the two approaches 
see section VII. 



IV. SOLUTION OF THE DYNAMICAL 
EQUATIONS, LIGHT-SCATTERING SPECTRA 

Due to the formulation of the problem in terms of 
Laplace-transforms of correlation-functions, we only have 
to invert matrices in order to solve (III. 4). The presen- 
tation of the results proceeds in three steps. First we 
solve the equations for helicities v = ±2, then v = ±1, 
and finally the v = part. This allows to go from the 
most simple coupling, namely none at all, towards the 
coupling of all variables. For notational convenience we 
abbreviate the Laplace-transforms of the dynamical vari- 
ables by C v sr := £{(p us \ Pv,r)} where v is defined by the 
context the equation is given and mostly omitted at all, 
and C pp :=£{(p\p)}. 

The helicities v = ±2 only concern orientational vari- 
ables, 

Cf 2 (z) = S sn (5 nr z + T nr (z)) D-^(z). (IV. 1) 
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where we have defined 



(IV.2) 



Due to the additional coupling of orientation to the mass 
current for v = ±1 one will find more diverse dynamics 
in these variables. Defining a damping f](z) by 



fj(z) := r) s {z) + zp, n (z)(A') nr D r ^(z)p, m (z). 
the solution for the oricntational variables reads 



where an amplitude 



2 C\ l a m (z)a s {z)m 2 v 2 
1 *+&CP ±1 fj{z) 



a m (z) := -^jJ, n (z)Q nr D r ^(z) 



(IV.3) 



(IV.4) 



(IV.5) 



we stick to the model of [4, 6] where Se := ap + Q is 
the sum of fluctuations of the mass-density and contri- 
butions of the orientational dynamics. The general di- 
electric tensor is a 3 x 3 matrix, thus a nine-dimensional 
object. Its representation in terms of irreducible repre- 
sentations of the rotation group is lying in a direct sum 
of the invariant subspaces with j = 0, 1, 2. The molecules 
of the liquid have been assumed to be optically inactive. 
This yields a symmetric dielectric tensor, thus having the 
j = 1 component vanishing. As discussed in section II the 
mass-density fluctuations will account for the scalar part, 
whereas the orientational variables, defined in (II. 3), give 
the tensorial component. The exact derivation of the 
scattering-spectra in cartesian coordinates can be found 
in [14]. Using the transformation laws between cartesian 
and spherical representations [13] one obtains up to an 
overall multiplicative constant the following expressions. 
The polarised spectrum is given by 



has been introduced. The v = components exhibit the 
most structure, for there are also couplings to the fluctua- 
tions of the mass-density. Introducing another damping- 
kernel 



rVV 



-f° 
> z -o.o 
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(1.2 + + o^2,2 



(IV.10) 



K(z) := Vb (z) + C 2 f)(z). 
one finds the density-density auto-correlator 



2 2 



q 2 c 2 



q 2 c 2 _ z 2 _ JU zK ( z ) 
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(IV.6) 



(IV.7) 



And the cross-polarised (in the literature this is also 
noted as 'depolarised') spectra reads 



rVH 




\ T7 ) 4,2 \- (IV.ll) 



^ ) C 2 ' — 2 



The autocorrelation-function of the orientational vari- 
ables is 

C% m {z) = Cf 2 (z) + 2 C 2 a m (z)a s (z)zm 2 v 2 
SmV ' SmW q 2 c 2 - z 2 - ^zK(z) 

And finally the correlation of density and orientation 
reads 



£° P ,o( z ) = C q 2 a m (z) 



2 2 

m v 



q 2 c 2 _ z 2 _ JlL zK iz) 
* run ^ ' 



(IV.9) 



A. Light-Scattering Spectra 

In order to describe a light-scattering experiment a 
model for the dielectric-tensor has to be assumed. Here, 



The dielectric tensor has been identified as an object Ssl 
where j = 0, 2. The according Laplace-transforms of 
its autocorrelation-functions and the limit of the imagi- 
nary part to the real axis, which is the actually measured 
quantity, have been abbreviated 



Cl k {u) := hm Im£{(4(t) | ej(0)>}(w + is), 



where ui € M and e > 0. The scattering-angle is denoted 
t9. The correlators of the oricntational components are 
obtained by multiplying those of the oricntational den- 
sities p^ m by the local polarisability tensor Q B n which 
without loss of generality we can assume to be diago- 
nal. Plugging in the expressions obtained in the previous 
section one finds the following 



tVV 



— lim Im 

3 e\0 



(-f-OB.flD-l.QB + ^I 

I lOz z 



Zq2 c 2_ z 2_^ zK{z) 



; C a-Q B 
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(IV.12) 



jVH _ 



- lim Im Q B { -SID 1 - cos 2 f -) q 2 C 2 



2 2 

m v 



z+^C 2 ±1 fi(z) 



a[z)a{z) W 



(IV.13) 
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NB: Bold symbols are the coordinate- free representation 
of the respective object, eq. Q ~ Q^. Both spectra ex- 
hibit a purely rotational background whose amplitude is 
independent of q. The second term of each spectrum is 
proportional to q 2 as expected by symmetry considera- 
tions [12]. In (IV. 12) this term represents a longitudinal 
phonon with velocity c, damped by the kernel K (z) and 
having a frequency-dependent amplitude which is also de- 
termined by the translation-rotation coupling hidden in 
a. The cross-polarised spectrum on the contrary shows a 
primarily diffusive shear mode with the viscosity C^ry. 
In a backscattcring geometry where $ = 7r this term van- 
ishes and the spectrum becomes purely rotational. Ef- 
fects of structural relaxation in the oricntational degrees 
of freedom are therefore accessible best in such a config- 
uration. If one assumes a characteristic behaviour of the 
memory-kernels at high and low temperature, one can 
derive limiting expressions for the spectra and the occur- 
ring phenomena. At high temperatures the frequency- 
dependence of the transport-coefficients can be neglected 
and the q 2 mode in I VH therefore becomes a purely dif- 
fusive shear mode, which is subtracted off the rotational 
background. This is the so-called Rytov-dip [18]. At very 
low temperatures one can model the structural relaxation 
of the system by non-ergodic correlation-functions [19]. 
The Laplace-transforms of such functions are negatively 



proportional to -. In this case the shear-mode becomes 
a propagating shear wave, showing the onset of solidifi- 
cation of the system near the glass-transition. Thus, the 
theory developed so far and the resulting expressions for 
the light-scattering spectra (IV. 12, IV. 13) are able, along 
with reasonably justified approximations of the memory- 
kernels, to explain the cross-over in the cross-polarised 
spectrum from one exhibiting the Rytov-Dip to the low- 
temperature case, where propagating shear-waves are ob- 
served [20]. 

If one neglects that there is a contribution of the 
translation-rotation coupling to the longitudinal phonon 
dissipation kernel K(z), the isotropic part of the po- 
larised spectrum, i.e., that part that is independent of the 
anisotropy of the molecules, is the longitudinal phonon 
with amplitude a, 

I iso :=aq 2 — hm Im - — . (IV. 14) 

3 s\o z „2 C 2 _ z 2 _ q_ zK r z \ 

^ van v ' 

In contrast to results of previous approaches, cf. [14], 
where I vv = Ii so + ^I VH ('& = n), we find an additional 
contribution to I vv due to the translation-rotation cou- 
pling, i.e., 



T vv 



I iso + -I VH ($ = ir)+q 



2 C 2 
— = lim Im< 

3V5 e\.o 



q 2 c 2 — z 2 — — zK(z) 

I 



—a 



V3 



a Q 



(2),S 



+ ■ 



(a-Q 



(2),BV 



>. (IV.15) 



The difference I vv - (\l VH {d = tt) + J iso ) has recently 
been measured experimentally and called the 'VV-Dip' 
[20] . This result shows that an approach without consid- 
ering the rotational dynamics or the translation-rotation- 
coupling is too simple and contradicted by experiments. 



V. POSITIVITY OF THE SPECTRA 

An application of the theory presented in the previ- 
ous sections to experiments is only possible when the 
memory-kernels are approximated by phenomenological 
functions. One problem, however, of any approximation- 
scheme is to find the admissible set of kernels that will 
lead to physical results or be able to identify those ap- 
proximations which yield unphysical results. For instance 
models for the memory-kernels should not lead to neg- 
ative scattering-spectra which renders positivity a use- 
ful criterion for the admissibility of an approximation- 
scheme. In this paragraph, we will discuss positivity of 
the correlation-matrix which solves (IH-4) and thereby 
prove that the spectra will be positive as long as the 
memory-kernels are positive semidefinite matrices. In 



section III, we have formulated the equations of motion 
as matrix equations for the Laplace-transforms of the re- 
spective correlation-functions. These equations are quite 
simple, they all have the form 

A-B = -ld, (VI) 

with A and B defined appropriately. Now let A be the 
space C" xn with positive and finite integer n, the space of 
all complex n x n matrices. With the terms of real and 
imaginary part of a matrix C we denote the following 
hermitian matrices, cf. [21], 

RcC := l{C + &) (V.2) 
ImC := ^(C-Ct), (V.3) 

where the dagger t denotes transposition and complex 
conjugation. A hermitian matrix A in A is said to be 
positive, i.e., A > 0, if its spectrum, which is real, also 
is nonnegative. Equivalently we can say it is positive 
semidefinite. Now let A, B e A and A - B = —Id. Mul- 
tiplying equation (V.l) with B^ from the left we find 
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B^-A-B = -B f and the adjoint equation B^- A^-B = -B. 
Now, subtracting these two one obtains 



B f ImA- B = lm.B. 



(V.4) 



It can be shown that for a positive element A the trans- 
formed element B^ ■ A ■ B is still positive [21]. Thus we 
know from the positivity of Imi? that Im^4 must be a 
positive element of A. In the following, we will prove 
that the positivity of the imaginary part of a correlation 
matrix C(z) is actually equivalent to the positivity of the 
imaginary part of the memory-kernel matrix (A | R'A) in 
the above sense. 



A. Application to the Mori-Zwanzig formalism 

According to equations (III. 4, V.4) we have to show 
that SMS > 0. All our considerations comprise a fi- 
nite dimensional subspace of the state-space. Without 
loss of generality we can therefore use orthogonal vari- 
ables, i.e., (Ai \Aj) — Sij (Ai | Ai) If the variables were 
not orthogonal we would have had to use a similarity 
transformation to make them orthogonal. However, the 
imaginary part of any matrix is hermitian, and hence the 
similarity transformation that diagonalises the imaginary 
part is unitary and will not influence any positivity prop- 
erties. In the orthogonal system, the matrix S^ is diago- 
nal and, because this holds for autocorrelation-functions 
in general, the diagonal elements are real and positive. 
The same is true for its inverse (Sij) -1 . The multiplica- 
tion by 5 -1 is therefore also irrelevant for the discussion 
of positivity. We now analyse the three different addends 
in equation (III. 4) and examine their imaginary parts. 

• Sij is real and symmetric. Thus ImSz = Slmz > 
or vanishes if we write z = ui + is and let e \ 0. 

• The imaginary part of the reactive contribution 
— i (Ai | Aj) is identically zero, because (Ai \ Aj) = 
-{Ai\Aj). 

• The last term represents the matrix of memory- 
kernels which is at our disposition. 

Collecting these properties and using the result about 
positive elements in C™ x " we have proved that with 
the memory- kernels having positive imaginary part, the 
correlation-functions will also have positive imaginary 
part. The converse also holds due to a theorem in [22]. 



B. Conditions on the Memory-Kernels 



memory- kernel matrix has positive imaginary part. 
We always assume the individual memory-kernels to 
be proper correlation- functions, i.e., be positive and 
bounded in the sense of Appendix C. Being given the 
matrices (III.7TII.9) we can extract their imaginary part, 
i.e., the imaginary part of the memory- kernel matrix 
(A | R'A). We will treat v = in detail here, the re- 
mainder works out by analogy, 



(A\R'A) V=0 = 



Its imaginary part is the matrix 



/0 





o\ 


f 


-iqCofin 





iqCoff 


q 2 i) 





\0 





0/ 



(V.5) 









Im f„ r 
iqC Im jl^ 
\0 







-iqC \m fi n 
q 2 Im?j 









(V.6) 



Since the matrix T nr is symmetric, the imaginary part 
can be taken element-wise. For the vector /2„ the same 
applies. Take a vector (a, b) where a e <C 2v+1 and b e C. 
Hence, positivity of (V.6) is equivalent to the inequality 



< (a b)" 



-iqC lmfi n 

r.2 



Im T nr 
K iqC Im q 2 Im fj 

= a* Im fa + b* (q 2 Im fj)b + 2 Re [b(iqC Im /2)a] . 

(V.7) 



This relation holds true for all b G C and in particular 
for the minimising value b min . Due to the convexity of 
b*(q 2 Im?7)6 in 6, the global minimum is the only critical 
value of the right-hand-side. Differentiating the above 
inequality with respect to b and equating to zero yields 



7 * 

min 



iqCo Im(/i)a 
q 2 lmf] 



Inserting into (V.7) we find that for positivity of the op- 
erator in question the memory-kernels need to fulfil 



Imr7(a* Imfa) - | Im(/2) ■ a\ 2 > 0, 



(V.8) 



which is both necessary and sufficient for the assertion 
to hold true. Condition (V.8) resembles the one found in 
[6]. However, this result as the application of the positiv- 
ity properties of Mori-Zwanzig-projections is much more 
general as it is derived using a method that does not rely 
on the specific properties of the system in question. 



We shall now apply the general result of the previ- 
ous paragraph to the equations of motion of the rota- 
tional dynamics in a super-cooled liquid. We will thus 
find necessary and sufficient conditions for the individ- 
ual memory-kernels of the theory such that the resulting 



VI. ENERGY FLUCTUATIONS 

In the preceding discussion, we have neglected temper- 
ature, i.e., energy-density, as a hydrodynamic variable. A 
complete theory of super-cooled molecular liquids should 
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include the hydrodynamic mode which is due to the con- 
servation of energy. In the available literature, fluctu- 
ations of the energy-density are usually left out of the 
discussion, mainly because due to the closeness of the 
Landau-Plazcek-ratio to unity the Rayleigh-line is sup- 
pressed if compared to the longitudinal phonon peaks. 
However, technological progress in measurement allows 
to detect the Rayleigh-line in light-scattering experi- 
ments [25, 26]. In order to provide a proper theoretical 
setting for these experiments, we briefly discuss how the 
theory can be extended to include energy-fluctuations. 
The very formal presentation of the projection-operator 
formalism allows an easy and transparent modification of 
the equations. 



conservation of energy, 

- ic°<d(q) = c° v Ce(q) = CQ p E K 

= q .jE_(ElA q .j- C QE p . (VI.3) 

(p\p) 

This quantity docs not define eigenvectors of the 
Liouville-operator for the eigenvalue q = as the mass- 
current and the momentum-current do. We therefore 
must relax the assumption, that memory-kernels can 
be approximated by their q — value. Indeed, we 
will actually find that the long-wavelength expansion of 
(0 | R'(z)Q) generates the heat-mode in the density-auto- 
correlator. 



A. Definition and Basic Properties 



B. Correlation- Functions involving Temperature 



In this section, we will introduce the notion of a ki- 
netic temperature whose fluctuations can be regarded as 
the variable being responsible for the Rayleigh-line in the 
spectra. According to [8] that part of the kinetic energy 
fluctuations which is orthogonal to the mass-density fluc- 
tuation is a suitable quantity. Most of the definitions and 
formal properties can also be found in [8] . Notice, that we 
decompose the total energy into E = E K + E p , i.e., ki- 
netic and potential energy fluctuations respectively, and 
label the kinetic specific heat at constant volume by c v . 
We define the temperature by 



c"e := Q P E 



K 



E 



K 



(E K \p) 
(P\P) ' 



(VLl) 



This kinetic temperature is normalised to 
c° v (0(q) | 6(g)) = k B T 2 . The total energy is con- 
served and thus we have an additional local equation of 
continuity, 



d t E = iq- J L 



(VI.2) 



Definition (VI. 1) indicates, that the temperature is a 
j = object, i.e., it behaves like a scalar variable un- 
der rotation. Exploring the possible static correlation- 
functions we find, that within the chosen set of hydro- 
dynamic variables in the long-wavelength limit tempera- 
ture will only couple to the mass-density and the poten- 
tial energy. Considering the microscopic form of these 
correlation-functions and the fact that the potential en- 
ergy only depends on positions, but not on momenta, re- 
veals that temperature is orthogonal to all variables that 
depend on positions only, in particular mass-density and 
potential energy. The projection-operator (III. 1) is there- 
fore trivially extended by adding the normalised kinetic 
temperature. For static correlators one may also replace 
Q p E p = QE P , where Q is a projection orthogonal to 
the extended hydrodynamic subspace, i.e., QO = 0. The 
time-derivative of the kinetic temperature is, using the 



Due to the decoupling of the equations of motion 
(III. 4) with respect to helicities, the only appearance 
of the temperature will be in the v = equations. 
By time-inversion parity the static correlation-functions 
(p | 6) = - (6 | p) = as well as (6 | 6) = 0. Using the 
equation of continuity for the mass-current (II. 8) 



(JO\Q) = J(P\Q P E K )+O(q 2 ). 



(VI.4) 



The correlation-functions involving the reduced 
evolution-operator R! and the thus obtained transport- 
coefficients are computed in the following. Beginning 
with possible correlations to the mass-density we find 

(p | R'Q) = -iq (Jo | R'Q) = ~iq (Jo \ QR'Q) = 0. 

Thus there are no correlations between mass-density and 
temperature. Next, we consider the correlations to the 
mass-current, 

(j \R'Q) = -iq(p\R'Q)+O(q 3 ) 

= ^-(p\R'£QE p )+0(q 2 ) 

= ^[z(p\R'E p ) + (p\Q {p) E p )]. 

The p in brackets in the last line may be put or not 
due to the possible replacement Q p E p = QE P . Using 
equations (III. 4, VI.4) we have a matrix element (cf. [12]) 



q ; [\P I 
qmk B T 2 



: (Jo | 6) + (Jo I R'Q) = \ [<p | Q P E) +z(p\ R'E P }} 

Co, 



(VI.5) 



which defines a generalised dynamical tension coefficient 
(3(z). Under the reduced time evolution (6 | R'Q) is an 
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autocorrelation-function. According to (A. 5) its expan- 
sion in powers of q reads 

(6 | R'Q) = zc v (z) + q 2 X(z) + 0(q 4 ). (VI.6) 

The generalised specific heat cy(z) is defined cy(z) := 
Cy + c v (z). The q 2 terms of the reduced temperature- 



autocorrelation-function define the generalised thermal 
conductivity A(z). For the definitions and terminology 
of these coefficients see also [8] . 

We can now compose the matrix M${z) and have 
therefore extended (III. 9) to a model which incorporates 
energy-fluctuations. Temperature is simply added in the 
fifth row and column of the correlation matrix. 



i/=0 



f zSnr 







iqCfx^ 








-iqCjin 
q 2 T] + zm 2 v 2 






2 2 

qm v 



2 2 

qm v 

qmk B T< 
3 






qmk B T 2 r> 

o 



(VI.7) 



(q 2 X + zc v ) J 



The solutions of the such modified equations of motion 
will only slightly differ from those which we have ob- 
tained in section IV. The autocorrelators of the orienta- 
tional density and that of the mass-density are modified 
in the longitudinal-phonon-dissipation-kernel K(z). One 
therefore has to replace 



1 v( \ 1 vt \ rnTf3{z) 2 

K(z) h-> K{z) - — , 

ran mn q z \{z) + zcy(z) 



(VI.8) 



introducing an additional wave-vector dependence in 
the damping K(z), which will therefore be denoted by 
the symbol K(z,q). The temperature-autocorrelation- 
function is 



q 2 X(z) + zc v {z) 



zq 2 T 2 m 2 v 2 (3(z) 2 



(q 2 \(z) + zc v (z)) 2 ( 



c 2 q 2 



^K(z,q)) ' 



(VI.9) 



Here, K (z, q) is the modified damping (VI.8) with the 
modification (VI.8). In the remainder of this section, we 
will discuss the question whether the modified phonon- 
propagator [c 2 q 2 — z 2 — ZjL-K(z, q)]^ 1 exhibits an addi- 
tional hydrodynamic mode, which is visible in the light- 
scattering spectra. We are not able to answer this ques- 
tion for the general solutions with the full frequency- 
dependence of the memory-functions. However, in 
the limit of high temperatures we can use a Markov- 
approximation for the retarded correlators. The modi- 
fication of the phonon-dissipation-kernel K(z, q) only af- 
fects the polarised spectrum I vv , cf. equation (IV. 12). 
We can therefore consider the autocorrelation of mass- 
density (IV. 7) which enters the expression for I vv using 
the modified phonon-dissipation-kernel K(z, q) to find 



the heat-mode, 

C{{p{q)\p{q,t))}{z) 



z — 



This correlation-function has the same functional form as 
the result for the density-auto-correlator that has been 
obtained in [8] and [12]. Using Markov-approximations 
in the high-temperature limit 

K(z,q)^ir}°, \(z)^i\°, 
c v {z) -> c v + izc v , (i{z) —> /3 + iz/3", 

with all constants being real and the corresponding limit 
of the memory-kernel matrix being positive definite, 
these authors have found three hydrodynamic modes, two 
sound waves and one heat-diffusion mode with poles at 



^sound — ^-C s q 



ia. 2 ^ + 0(q% 



zhoat = -iq D T . 



The approximation for small q is due to [23]. The adi- 
abatic sound velocity is defined as c 2 := c 2 + raT(3 2 /cv 
and the dissipation-constants are given by 



D 



rf 
mn 
c 



mT(i 2 
c v 



c 2 c cv f3 



2 A 



c 2 c v ' 



,(VL10) 
(VI.11) 



In contrast to the sole viscous and thermal damping of 
the sound waves by the first and the second term in 
(VI. 10) that is known from hydrodynamics [16], T; has 
two additional contributions. These stem from the use 
of a kinetic temperature in our formulation. Since the 
total energy is conserved, there are contributions in the 
correlation-functions showing how kinetic energy is trans- 
formed into potential energy This effect is irreversible 
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and therefore displayed by two additional terms in the III. For v = for instance it has the form 



damping of the sound-waves iy The term with c v de- 
scribes the irreversible transfer of kinetic to potential en- 
ergy at constant volume. Heating of the liquid due to 
mechanical work, i.e., the transformation of clastic en- 
ergy into potential energy, is accounted for by the term 
with j3". 



VII. COMPARING TO PREVIOUS RESULTS 

Recently, a number of theoretical results for the prob- 
lem of orientational dynamics and its application to light- 
scattering have been published [4-6, 12]. The theoretical 
description given in these papers is comparable to the 
approach pursued in this work. However, the remark- 
able difference lies in the fact, that they directly use the 
decomposition of the dielectric tensor into its irreducible 
components effl as basis-vectors for the hydrodynamic 
subspace, whereas we have preferred to project the dy- 
namics onto the so-called tensorial densities p^ m which 
by definition are closer related to the rotational dynamics 
of the molecules; see equations (II. 3, II. 4) for the relation 
of the two quantities. The purpose of this section there- 
fore is to relate our theory to the existing results by first 
presenting a theory, in the following called the scalar the- 
ory, which projects onto the variables Q M and then mod- 
ify the equations of section III to yield exactly the same 
set of equations. It is obvious that the main concern is 
the modification of the memory-kernels and the question 
whether the newly derived kernels can still be interpreted 
as memory- kernels, i.e., as correlation- functions of a cer- 
tain stochastic process. 

For the discussion of the latter we will define the no- 
tion of a correlation-function in an analytically tractable 
way as the characteristic function of a certain probabil- 
ity measure. Then properties of the Laplace-transform of 
such objects are used to prove that the functions which 
will appear in place of the memory-kernels of the scalar 
theory can indeed be interpreted as memory- kernels. The 
mathematical theory which we will have to employ is dis- 
cussed in [22, 24]. A brief overview regarding terminology 
can be found in Appendix C. 

If one directly considers the irreducible components of 
the dielectric tensor as part of the basis of the hydrody- 
namic subspace in an application of the Mori-Zwanzig- 
technique, a suitable projection-operator will have the 
following form 



V(q) :- 



IQmXQmI , \Q»)(QA , , \p)(p\ 



+ 



n 



+ 



+ 



(Jo\Jo) (p\p)' 
(VII.l) 

where e*- -* <~ p, effl <~ Q^, and proper summation over 
\x is implied. The normalisation-factors are defined in a 
similar fashion as in section II. Accordingly, the matrices 
M v are derived by calculations following those of section 
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-iqCofb 
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qm v 
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(VII.2) 

where the definition of the memory-kernels formally is 
the same as before. 

In the following discussion we will begin with giving 
results for the general case of dimensional reduction in 
Mori-Zwanzig-applications and then specialising to the 
case of orientational dynamics where appropriate. 



A. Some General Results 

Consider the equation for a matrix of correlation- 
functions as given by Mori and Zwanzig in the form [16] 



C(z) = -S(z-lj + M(z)Y 



(VII. 3) 



C(z) represents the correlation-matrix, S the normalisa- 
tions, iv the libration frequencies, and M the necessary 
memory-kernels. Without loss of generality we can as- 
sume that M(z) is of order less than i. The general 
representation of analytic functions with positive imagi- 
nary part, as M is such an object, has a term linear in 
z, one being constant and the remainder being of order 
less than i, see [22]. The linear and constant parts can 
be absorbed into S and u>. As can readily be seen equa- 
tion (VII. 3) is equivalent to (III. 4). The reduction to the 
scalar theory is done by multiplication with suitable ma- 
trices a from the left and the right. These produce the 
associated quadratic form for each block-matrix. The 
reduced correlation-matrix obeys the equation 

C Icd (z) = a T C(z)a = -a T S (z - uj + M(z)) -1 a, 

(VII.4) 

which can be cast into 

C Yc d{z) = -5 rc d (z - Ld Icd + M Icd {z)) 1 . 

thus defining reduced quantities S re d, w roc i, and M rod . As 
before we can assume M rc( j to be of order less than - , be- 
cause according to our discussion of the positivity prop- 
erties of correlation-matrices, M roc j is an analytic func- 
tion with positive imaginary part. We can use the high- 
frequency expansion of the correlator to gain more infor- 
mation. The zeroth moment of the representing measure 
of M(z), cf. Appendix C, will be denoted mo. 



C re d(z) ~ a 1 



z z 



Suj S / 9 \ 



(VII.5) 

The high-frequency expansion of C ro d with the reduced 
quantities is 



C re d(z) 



— SVod SrcdO^rcd 5* ro d 
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Comparing the coefficients we find 



where u> = S 1 £l is used. Thus 



Sred 


= a T Sa 


(VII.6) 




^rcd 


= a T Sua 


(VII.7) 


n = Slu 


W re d 


= ^rod^rad 


(VII.8) 
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(mv) 2 q 



and obtain an equation for the masses of M ved . 

m 0,rod = S'rcd [(Sm )rcd ~ (^)rcd + ^rcd^rcd] 



(VII.9) 



where (S , m ) rc d = a^moa and (£lcj) ro d = a^flua. 
Thus, the measure which is representing the correlation- 
matrix M ro( j has finite mass and renders the matrix a 
correlation- function in the sense of Appendix C. We 
note by the way, that the definitions of S re d and w ro d are 
in accordance with the terminology and notation of [6]. 
That is normalisation and libration frequencies survive 
the reduction without any loss. The only changes will 
appear in the memory- kernels. From a physical point 
of view this is quite obvious. Each of the rotational 
variables has its own libration frequency. Contracting 
these generates beats. The frequency u> ro d is the carrier, 
the remaining oscillations are put into a friction term, 
i.e., memory- kernel. The fact that this extension of the 
memory-functions does not alter their property of being 
the correlator of some stochastic process is shown by the 
above calculation. 



B. Number of Memory- kernels 

The remaining part to prove is that the number of 
memory-kernels in the reduced theory equals the number 
of memory-kernel matrices of the vector-valued theory. 
We do not bother with treating all helicities separately, 
but give the results for the most difficult case, v = 0. 
Then the rest is obvious. Consider the memory-kernel 
matrix of the vector- valued theory, 



M(z) = 
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To compute the number and position of the non- 
vanishing elements of the reduced memory-kernel matrix 
we start by explicitly computing the respective masses 
™o,rod, defined in cq. (VII.9). The matrices involved in 
its computation are the following. 



S = 



and 
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and hence 



Qui = 







(mv) 2 (, 







-(mvq) 2 

-(mvq) 2 ■ 



The matrix fi re dW r ed has a form analogous to this. Hence 
we infer that the only non-vanishing element of the dif- 
ference fircd^rod — (f^)rcd is the (2, 2)-element. S~ d 
is diagonal, and will not change this fact. Since we 
know the form of mo, we find that mo >re d has only the 
elements (2,2), (2,3), (3,2), and (3,3) non- vanishing. 
It is a fact from function-theory, that the Stieltjes- 
transform of a measure has zeros on the diagonal, when- 
ever the mass- matrix vanishes there [27]. The imag- 
inary part of the memory-kernel matrix M rc d is posi- 
tive definite. For positive definite matrices A we know 
\A jk \ < A^A 1 ^. And therefore (Im M rod )ij = for all 
(i,j) in {(l,fe), (k, 1), (4,fc),(fc,4) I 1 < k < 4}. However, 
the memory-kernels are analytic functions and thus the 
real part of these matrix-elements can only be a constant 
which by the assumptions on M re d is zero. Hence the 
reduced memory-kernel matrix has the structure 



M rod = 



/0 0\ 

A B 

C D 

V0 0/ 



This means we have four memory-kernels at most! Now, 
we have to prove, that it is actually three distinct 
memory-kernels, in particular we show that B is pro- 
portional to C. We can give an explicit expression for 
M rcd (also cf. [22] pp. 111.), 

M rcd (z) = -z + w rcd + S red C~l(z). 

The first term only affects the diagonals A and D and 
the matrix w ro d has no nonzero elements in the positions 
of interest. So we only have to care about the inverse of 
the reduced correlation matrix. The key of the proof are 
the specific symmetry properties of C~^. The reduced 
correlation matrix has the following form 



red 



Cs 

c c ■ c c 

t^u.s *"i/,s t^u.v ^v.p 

^p,s ^p.s ^P-V ^P^pf 



c ■ c 
c ■ ■ c ■ 



s,p 



The subscripts of C denote the various variables which 
are correlated in order to form the respective correlation- 
function. The subscript s denotes the rotational degrees 
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of freedom, s the associated currents, v the mass-current, 
and p the mass-density. We remember that for irre- 
ducible representations of the rotation group the inter- 
change of variables in correlation-functions obeys equa- 
tion (A. 7). Thus the parities under exchange of variables 
form the following matrix 




Now, how do the symmetries transfer to particular ele- 
ments of the inverse matrix? In general, the inverse A^ 1 
of a matrix A has elements 



^4-1 = _L_ 

i k dct A 



kj- 



With Akj we denoting the determinant of the matrix 
which is identical to A except that the k th row and the 
j th column are missing. In particular, we are interested 



in the elements A } and An}. Thus A 2 



-A, 2 . This 



i 23 ai,u ^132 . n, uo si 2 3 — -^32- 

shows that the elements of interest in the inverse of the 
reduced correlation matrix are related to each other by 
real scalars. In more detail 



M lcd (z) 23 = -^—C rcd {z) 23 

S ^rcd 



(VII. 10) 



M Icd (z) 32 = 



-1 



kBTm 



n ied M ied ( z ) 23 . (v11.1i) 



This shows the formal equivalence of the vector-valued 
theory and the scalar theory. To proceed from the first 
to the second one only has to replace the block-matrices 
by scalars. The memory-kernels retain their property of 
being correlation-functions. And due to positivity even 
the necessary condition on the memory-kernels [6] turns 
out to be the same. 



the complete rotational dynamics by using tensorial den- 
sities derived from Wigner-functions. The such obtained 
equations of motion can easily be solved using Laplace- 
transform methods and are used to obtain expressions 
for the light-scattering spectra both in the polarised and 
the depolarised configuration. We were able to repro- 
duce the main features of the experimentally measured 
results, i.e., the Rytov-Dip at high temperatures and 
its vanishing in favour of propagating shear-waves af- 
ter cooling. The importance of the translation-rotation 
coupling not only for the Rytov-Dip has been demon- 
strated by showing that the intensity in the VV config- 
uration is not I times that of the VH configuration in 
backscattering-geometry plus the isotropic contribution, 
thus contradicting the classical result for simple liquids 
[14]. This difference, the so-called W-dip, due to the 
translation-rotation coupling has recently been measured 
[20] . The role of energy-fluctuations for the spectra has 
also been discussed. The theoretical model has been 
extended by also considering hydrodynamic modes of 
energy-fluctuations by similar techniques as used before. 
The appearance of an additional line in the scattering- 
spectrum due to heat, the Raylcigh-line, has also been 
shown. Finally we were able to derive equations of mo- 
tion from our more general theory which are formally 
equivalent to those obtained by only considering the ir- 
reducible components of the dielectric tensor itself in the 
projection-operator formalism. Using results from the 
theory of positive functions and the classical moment 
problem of probability theory we were able to charac- 
terise the Laplace-transforms of correlation-functions as 
those positive functions which satisfy a certain growth 
condition. This result enabled us to identify the func- 
tions which appeared in place of the memory-kernels in 
the simpler theories as correlation-functions, hence sup- 
porting their interpretation as memory-kernels as well. 



VIII. CONCLUSIONS 

In the present work, we have given a complete and com- 
prehensive account of the long-wavelength rotational dy- 
namics of a super-cooled molecular liquid and computed 
expressions for the light-scattering spectra of such sys- 
tems. The theory is based on the concept of generalised 
hydrodynamics which views the transport-coefficients as 
frequency-dependent quantities, and an application of 
the projection-operator-tcchniquc of Mori and Zwanzig. 
The frequency-dependence of the transport-coefficients 
which appear as integral-kernels in the generalised 
Langevin equations allows both the description of the 
liquid in the high-tempcraturc regime, where such de- 
pendence is negligible and the physical behaviour is that 
of a simple liquid, and the modelling of the solidifica- 
tion of the liquid near the glass-transition as it is known 
from mode-coupling theory. In contrast to previous ap- 
proaches that only consider a decomposition of the di- 
electric tensor into its irreducible components, we model 



APPENDIX A: SELECTION RULES 

Due to the rotational invariance of the Hamilto- 
nian function we can state a couple of selection rules. 
These will greatly simplify the evaluation of correlation- 
functions. We always assume that q is given parallel to 
the z-axis, i.e., q — qe$. Most of these rules can be de- 
rived easily, for a more complete presentation see [10]. 
Given two variables A and B lying in irreducible repre- 
sentations of rank i and j respectively it follows 

(A^(q)\B^(q))=S^(q). (A.l) 

Using rotation-matrices and playing around a little one 
finds 

S^(qT = (-iy +3 S^(q) (A.2) 

and 

s?( g ) = (-i)'M4n(?). ( A - 3 ) 



where e q denotes the parity under spatial reflection for 
each variable. In the special case q = an even stronger 
selection rule applies 

(Af\B<j))=5^5 ij S i . (A.4) 

If we expand the correlator (A.l) in powers of q the first 
non- vanishing term will be of order \i — j\. Two useful 
examples are 

S™(q) = S 00 + O(q 2 ) (A.5) 

and 

Sf{q)=0{q 2 ). (A.6) 

The Hamiltonian function is invariant with respect to 
time inversion. Hence we can assign time-inversion- 
parities e\ to the dynamical variables. And therefore 
we find for the behaviour of dynamical correlators under 
exchange of variables 

<4%t)|B<>'>(g)) = (-1)^6^6% (B^(q,t)\A^(q)). 

Time inversion parities and the tensorial rank of dynam- 
ical variables are preserved under time evolution. 

APPENDIX B: NORMALISATION 

In this Appendix we show the invcrtibility of the nor- 
malisation matrices S mn and fl mn as they have been de- 
fined in section II. Due to the simplicity of the proof 
we can deal with general j here. The case needed in the 
present paper is j = 2 only. The matrix S mn as defined 
in equation (II. 6) is a real and symmetric matrix on a 
space C 2j+1 by property (A. 2). Thus one may diago- 
nalise it. What is left to prove its invertibility is, that 
none of the eigenvalues vanishes. A sufficient condition 
is, that the associated linear operator is positive definite. 
Taking an arbitrary vector A e C 2:,+1 wc compute the 
scalar product 

^y m S mn Xn = (^2,KnP*m X nPn) = X mP m \ 2 ) > 0. 

ran van va 

Since the Kubo-scalar product is non-degenerate we only 
need to look whether there is a certain A to make the sum 
vanish. The proof goes indirectly: Let ^ A e C 2 - 7 " 1 " 1 
such that J2 m ^vaPm = 0. Recalling the definition of 
p^ m (II. 4) this condition is equivalent to 

N 

£e^£©C»(Q a )A m = 0. 

a — 1 va 

We assume that all the R a are different and since the 
exponentials are orthogonal we can project onto the con- 
tribution of one single molecule a, 

^D«(n a )A ra =0. 

m 
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Multiplying from the right with T>^ _ n (Ct a )* and inte- 
grating over the group we can exploit the orthogonality 
between the representations [28] and find A n = 0. Since 
n was arbitrary A = and the proof is given. 

To prove that the normalisation matrix of the orienta- 
tional current is invertible, is more complicated. Using 
angular momenta instead of Eulerian angles we formulate 
the orientational kinetic energy as 

a i=\ 1 

Some algebra yields 

n mn = ^±l fcB T^ J dn©o>(n)Mx>W(n). (B .2) 

k 

This is nothing but the quantum mechanical asymmetric 
rigid rotor whose solution can be taken from [29]. Like 
the matrix S mn , the normalisation-matrix of the orien- 
tational currents Vl mn is real and symmetric. Hence it 
is diagonalisablc. We only need to show its positive def- 
initeness. All involved operators appear quadratically 
and because the Li are self-adjoint for i = x,y, z, their 
eigenvalues are real and the eigenvalues of their squares 
are nonnegative. Consider that contribution to the ma- 
trix which is proportional to L 2 . This operator is already 
diagonal and its matrix-elements are 

(L 2 ) mm oc m 2 with m e n Z. 

Let A m <G C 2j+1 be the eigenvectors of L z . The opera- 
tor L 2 has a non- vanishing kernel, namely K :— lin{Ao}. 
Obviously, using the casimir operator L 2 , the identity 
L 2 - L\ = Ll + L 2 y holds and \* {L 2 X + L 2 ) 00 A = 
\* (L 2 ) 00 Xo > 0. When K is the kernel of Vl mn , we can 
either have Aq(L 2 )ooAo = or the respective quadratic 
form of L 2 vanishing, but not both, because their sum 
is positive. The only chance to find a zero eigenvalue 
would have been to have a common kernel of all three 
operators. The above relation shows that this cannot be 
accomplished and hence O m „ is positively definite and 
therefore invertible. 



APPENDIX C: POSITIVE FUNCTIONS 

To accomplish the reduction of the presented theory of 
light-scattering in super-cooled liquids to a scalar theory 
comparable to [4, 6], some mathematical definitions and 
results of the theory of positive functions are needed. 
We define a correlation-function to be the characteristic 
function of a probability measure [30-32]. 

0(t) := J e- M d^{\), (C.l) 
n 
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where without loss of generality J dfi(X) = 1. One can 
easily compute its Laplace-transform 

= / x^-z MX) - (C ' 2) 

This means that the Laplace-transform of a correlation- 
function is the Stieltjes-transform of a probability mea- 
sure. In particular, it is an analytic function mapping the 
complex upper half-plane into itself. In connection with 
the solution of the classical moment problem [22] , Ham- 
burger and Nevanlinna have proved that for an analytic 
function / mapping the upper half-plane into itself being 
the Stieltjes-transform of a finite measure it is necessary 
and sufficient that 

sw\yf(iy)\ < oo. (C.3) 

y>l 



The moments Sk of the measure jj,(X) are given by the 
coefficients of the large z expansion of /, 



So Si s 2 

z z 2 z 3 



Often the moments of the representing measures are 
called the moments of /. Thus an analytic function 
mapping the upper half plane into itself is regarded the 
Laplace-transform of a correlation-function if and only if 
the above criterion (C.3) on the growth of / holds true. 
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